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We compare variable number schemes for heavy quark electroproduction. 



Heavy quark production has been a major topic of 
investigation at hadron-hadron, electron-proton and 
electron-positron colliders. Here a review is given of 
some topics which are of interest primarily for electron- 
proton colliders. We concentrate on this reaction be- 
cause a theoretical treatment can be based on the 
operator product expansion, and also because data 
are available for deep-inelastic charm production at 
HERA. How all this relates to Fermilab experiments 
will be discussed at the end. 

In QCD perturbation theory one needs to introduce 
a renormalization scale and a mass factorization scale 
to perform calculations. We choose both equal to /i^, 
which will be a function of and the square of the 
mass of the charm quark w? . At small \j? ^ where 
kinematic effects due to quark masses are important, 
the best way to describe charm quark production is 
via heavy quark pair production where parton den- 
sities are used for the initial state light mass u,d,s 
quarks and the gluon and m only appears in the heavy 
quark coefficient functions (or partonic cross sections) 
like H^'^'^\z,Q'^ .nn^ , y?), etc., [ |l|. Here the super- 
scripts refer to their flavor decomposition and the or- 
der in perturbation theory, while the subscripts refer 
to the projection i — 2, L and the partonic initial state. 
The arguments refer to the partonic Bjorken variable 
z = Q"^ /{s + Q^) and to the fact that these functions 
depend upon invariants and scales. The renormal- 
ization necessary to calculate these NLO expressions 
follows the CWZ method [ |2j . The symbol H refers 
to those coefficient functions which arc derived from 
Fcynman diagrams where the virtual photon couples 
to a heavy quark line. Analytic expressions for these 
functions are not known. Numerical fits are available 
in [ H . Asymptotic expressions for » m? are avail- 
able in [ ^ which contain terms like \t?{Q^ /m^) and 
ln((5^/m^) ln((5^//i^) multipled by functions of z and 
demonstrate that these functions are singular in the 
limit that m ^ 0. There are other heavy quark coeffi- 
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cient functions such as L^^ (z,Q'^,'m^), which arise 
from tree diagrams where the virtual photon attaches 
to the an initial state light quark line, so the heavy- 
quark is pair produced via virtual gluons. Analytic 
expressions for these functions are known for all z, 
and m^, which, in the limit ^ contain powers 
of ln{Q^/m^) multiplied by functions of z. The three- 
flavor light mass MS parton densities can be defined 
in terms of matrix elements of operators and are now 
available in parton density sets. This is a fixed or- 
der perturbation theory (FOPT) description of heavy 
quark production with three-flavour parton densities. 
Notice that due to the work in [ |^ the perturbation 
series is now known up to second order. In regions 
of moderate scales and invariants this NLO descrip- 
tion is well defined and can be combined with a frag- 
mentation function to predict exclusive distributions 
[ 1^ for the outgoing charm meson, the anti-charm 
meson and the additional parton. This NLO mas- 
sive charm approach agrees well with the recent D- 
meson inclusive data in [ ^ and [ Q . Let us call the 
charm quark structure functions in this NLO descrip- 
tion Fl;^^^'^{x,Q^,m'^,nf = 3), where i = 2,L. 

A different description, which should be more ap- 
propriate for large scales, where terms in m? are neg- 
ligible, is to represent charm production by a parton 
density /c(a;,/x^), with a boundary condition that the 
density vanishes at small values of /i^ . Although at first 
sight these approaches are completely different they 
are actually intimately related. It was shown in [ 
that the large terms in \n{Q^ /m^) which arise when 
3> m^, can be resummed in all orders in pertur- 
bation theory. In this reference all the two-loop cor- 
rections to the matrix elements of massive quark and 
massless gluon operators in the operator product ex- 
pansion were calculated. These contain the same type 
of logarithms mentioned above multiplied by functions 
of z, (which is the last Feynman integration parame- 
ter). After operator renormalization and suitable re- 
organization of convolutions of the operator matrix el- 
ements (OME's) and the coefficient functions the ex- 



pressions for the infrared-safe charm quark structure 
functions Fi,c{x,Q'^ ,m'^ , A), become, after resumma- 
tion, convolutions of hght mass four-flavor parton co- 
efficient functions, commonly denoted by expressions 
like clf\Q^/fi'^) (available in [ §, [ with four- 
flavor light parton densities, which also include a charm 
quark density fc{x,^^). Since the corrections to the 
OME's contain terms in ln{Q^/m'^) and ln(m^/^^) as 
well as nonlogarithmic terms it is simplest to work in 
the MS scheme with the scale /i^ — rn^ for < m? 
and ^2 = ^2 _^ q2(^ „ mVQ2)2/2 for Q2 > ^2 ^^^^ 

discontinuous matching conditions on the flavor densi- 
ties at /i^ = m? . Then all the logarithmic terms vanish 
at Q2 = = 77j2 and the non- logarithmic terms in the 
OME's are absorbed into the boundary conditions on 
the charm density, the new four-flavor gluon density 
and the new light-flavor u,d,s densities. The latter are 
convolutions of the previous three-flavor densities with 
the OME's given in the Appendix of [ ||]. Hence we 
have a precise description through order of how, in 
the limit to — > 0, to reexpress the f'EXACT(^^g2^^2) 
written in terms of convolutions of heavy quark co- 
efficient functions with three-flavor light parton densi- 
ties into a description in terms of four-flavor light-mass 
parton coefflcient functions convoluted with four-flavor 
parton densities. This procedure leads to the so-called 
zero mass variable flavor number scheme (ZM-VFNS) 
for ^i.c(a;, Q^, A) where the m dependent logarithms 
are absorbed into the new four-flavor densities. To im- 
plement this scheme one has to be careful to use inclu- 
sive quantities which are coUinearly finite in the limit 
TO — > and A is an appropriate parameter which en- 
ables us to do this. In the expression for Fi^c there is a 
cancellation of terms in \t[^{Q'^ /m?) between the two- 
loop corrections to the light quark vertex function (the 
Sudakov form factor) and the convolution of the densi- 
ties with the soft part of the L2,g-coefficient function. 
This is the reason for the split of L^^g into soft and hard 
parts, via the introduction of a constant A. Details and 
analytic results for Lf ^^"^ and Lf^^^ are available in 
[ |ll| . All this analysis yielded and used the two-loop 
matching conditions on variable-flavor parton densities 
across flavor thresholds, which are special scales where 
one makes transitions from say a three-flavor massless 
parton scheme to a four-flavor massless parton scheme. 
The threshold is a choice of /j, which has nothing to 
do with the actual kinematical heavy flavour pair pro- 
duction threshold at Q'^{x~^ - 1) = 4m2. In [ D,[ 
H it was shown that the Ff^f^'^'^{x,Q'^,m'^,nf = 3) 
tend numerically to the known asymptotic results in 
Ftc™^{x,Q^,m^,nf = 3), when » to^, which 
also equal the ZM-VFNS results, which were called 
Ff^^{x,Q'^,nf = 4). The last description is good for 
large (asymptotic) scales and contains a charm density 



fc{x,fi^) which satisfies a specific boundary condition 
at /i^ = m^. 
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Figure 1. The charm quark structure functions 
Ff^^^'^ijif = 3) (solid line) F^f^[nf = 4), (dot- 
dashed line) F^f^^{nj = 4), ' (dashed line) and 
Fl^^ijij = 4), (dotted fine) in NNLO for x = 0.005 
plotted as functions of Q^. 

s most of the experimental data occur in the kine- 
matical regime between small and asymptotic 
a third approach has been introduced to describe the 
charm components of Fi{x,Q'^). This is called a vari- 
able flavor number scheme (VFNS). A first discussion 
was given in [ where a VFNS prescription called 
ACOT was given in lowest order only. A proof of fac- 
torization to all orders was recently given in [ Q for 
the total structure functions Fi{x,Q^), but the NLO 
expressions for Fi_c{x, , rn?) in this scheme were not 
provided. An NLO version of a VFNS scheme has been 
introduced in [ |ll| and will be called the CSN scheme. 
A different approach, also generalized to all orders, was 
given in [ ||,[ which is called the BMSN scheme. 
Finally another version of a VFNS was presented in [ 
p^ , which is called the TR scheme. The differences 
between the various schemes can be attributed to two 
ingredients entering the construction of a VFNS. The 
first one is the mass factorization procedure carried 
out before the large logarithms can be resummed. The 
second one is the matching condition imposed on the 
charm quark density, which has to vanish in the thresh- 
old region of the production process. All VFNS ap- 
proachs require two sets of parton densities. One set 
contains three-fiavor number densities whereas the sec- 
ond set contains four-flavor number densities. The sets 
have to satisfy the MS matching relations derived in [ 
p. Appropriate four-flavor densities have been con- 
structed in [ |l^ starting from the three-flavor LO and 
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NLO sets of parton densities recently published in 
f6|. 
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Figure 2. The charm quark structure functions 
j^EXACT(^^ = 3) (soUd hne) F^f^inf = 4), (dot- 
dashed line) Fl^^^{nf = 4), (dashed line) and 
Fl^^{nf = 4), (dotted Hne) in NNLO for x = 0.005 
plotted as functions of Q^. 



Since the formulae for the heavy quark structure 
functions are available in [ |l^ we only mention a few 
points here. The BMSN scheme avoids the introduc- 
tion of any new coefficient functions other than those 
above. Since the asymptotic limits for ;» of all 
the operator matrix elements and coefficient functions 
are known we define (here Q refers to the heavy charm 
quark) 



F, 



BMSN 



{x,Q^,m'^,A,nf = 4) 



F, 



EXACT 



{x,Q'^,m'^,A,nf = 3) 



ASYMP 



{x, 



,A,nf = 3) 



-^^;^F(x,g^m^A,n/ = 4). (i) 

The scheme for Fpg^ introduces a new heavy quark 
OME Aqq^^'(z,^Vto2) g^^^ coefhcient functions 

Hf^Q'''^\z,Q^/m^) because it requires an incom- 
ing heavy quark Q, which did not appear in the NLO 
corrections in [ The CSN coefficient fuctions are de- 
fined via the following equations. Up to second order 
we have 

Q2, 



C, 



A 



CSN, SOFT, NS, (2) 
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with the virtual term the second order Sudakov form 
factor. The other CSN coefficient functions are defined 
by equations like (we only give one of the longitudinal 
terms for illustration) 

,CSN,S,(l)/<3^ Q^" 



C 



L.g 



CSN,NS,(0) 



.S,(l)._Af_ 



CSN,NS 



'^i^). (3) 



with '"'''^^^ = 4m2/g2. The CSN and BMSN 
schemes are designed to have the following two prop- 
erties. First of all, suppressing unimportant labels, 
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(-/ = 4) = i^T^(n, = 4) 
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Since /q(i 
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^ (see [0]) this condition can be 



only satisfied when we truncate the perturbation series 
at the same order. The second requirement is that 



lim F^ 



BMSN 



("/ = 4) 



hm ^;CSN(„^_4) 



hm i^™F(n^.=4).(5) 



The only differences between the two schemes arises 
from terms in rri^ so they may not be equal just above 
— 777,2 ^j^jg turns out to be the case for the longi- 
tudinal structure function, which is more sensitive to 
mass effects. We plot in Fig.l NNLO results for the 
Q2 dependence of F|xacT(^^ ^ 3)^ i^2^f^("/ = 4), 
Fi^^^inf = 4), and ^^^^(n/ = 4) at x = 0.005. 
This figure shows that the results satisfy the require- 
ments in Eqs.(4) and (5). The ZM-VFNS description 
is poor at small Q^. In Fig. 2 we show the results for 

= 3), F^f (n/ = 4), FBMSN(„^ ^ 
and Fl^^{nf = 4) at a; = 0.005. We see that the 
CSN result is negative and therefore unphysicsl for 
2.5 < < 6 (GeV/c)2 which is due to the term in 
Am? IQ^ and the subtraction in Eq.(3). 

One way this research work is of relevance to Fermi- 
lab experiments is that it produces more precise ZM- 
VFNS parton densities. Such densities are used ex- 
tensively to predict cross sections at high energies, for 
example for single top quarks. Therefore the previ- 
ous work on four-flavor parton densities has been ex- 
tended in [ ^ to incorporate the two-loop discontin- 
uous matching conditions across the bottom flavour 
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threshold at /i = mb and provided a set of five- 
flavour densities, which contains a bottom quark den- 
sity fb{x,^'^). The differences between the five-flavor 
densities and those in [ ^ and [ ^ are also discussed. 
Results for deep-inelastic electroproduction of bottom 
quarks are presented in [ ^ . 
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